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Previous studies of lasers and nonlinear resonators have revealed that the polarisation degree 
of freedom allows for the formation of polarisation patterns and novel localized structures, such 
as vectorial defects. Type II optical parametric oscillators are characterised by the fact that the 
down-converted beams are emitted in orthogonal polarisations. In this paper we show the results of 
the study of pattern and defect formation and dynamics in a Type II degenerate optical parametric 
oscillator for which the pump field is not resonated in the cavity. We find that traveling waves 
are the predominant solutions and that the defects are vectorial dislocations which appear at the 
boundaries of the regions where traveling waves of different phase or wave-vector orientation are 
formed. A dislocation is defined by two topological charges, one associated with the phase and 
another with the wave-vector orientation. We also show how to stabilize a single defect in a realistic 
experimental situation. The effects of phase mismatch of nonlinear interaction are finally considered. 



I. INTRODUCTION 

Most studies of pattern and localized structure P-|l7|l formation in optical systems consider light with a fixed 
linear polarisation. However, the vectorial degree of freedom of light, associated with a space- and time-dependent 
polarisation, leads to a very interesting phenomenology, as previously shown for Kerr cavities [fl8|-p0[ , sodium vapors 
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]21| . lasers |2^-|2j| and intra-cavity second harmonic generation [25[j . This degree of freedom is also relevant from 
the point of view of information encoding and processing, because it yields a tool for the control of these structures. 
Optical vortices have been predicted to occur in optical systems such as lasers [ p6| and they have been experimentally 
observed in lasers |2^] and photo- refractive resonators |l2|,^8|,[29| . More recent theoretical studies revealed the existence 
of vectorial defects in optical systems . 

Among optical systems, optical parametric oscillators (OPO's) have been the object of intense theoretical study 
in recent years, which has revealed the existence of patterns pf]|- ffl[ and localized structures Jio| , p^|p^Jl7| , ^ -p5[ ; few 
examples of defect dynamics have also been published |l7]|58|,[l6| . Finally, experimental results on transverse pattern 
formation have been recently presented An OPO consists of a ring (or Fabry-Perot) resonator containing a 

nonlinear quadratic medium, which performs a parametric down-conversion of an injected laser beam (optical pump) 
at frequency u> p . Two new fields are generated in the crystal, the signal and the idler, at frequencies uj s , u>i respectively, 
such that lo p — oj s + u>i, the last relation indicating the conservation of photon energy. The nonlinear interaction also 
requires photon momentum to be conserved, i.e. the phase-matching condition k p — k s — ki — 0, which is usually 
obtained by exploiting polarized beams and crystal birefringence. This effect often implies that a transverse walk-off 
between the different polarisation components may occur. The consequences of walk-off have also been considered 
and convective instabilities and noise-sustained patterns found j3^J^,[w||49|,^0| . Most of these theoretical studies have 
been performed in the context of a scalar approximation, which is valid for the so called Type I OPO's; in the Type 
I interaction the signal and idler fields have the same polarisation, orthogonal to the pump beam. However, there 
are not many studies involving Type II OPO's where the polarisation degree of freedom is taken into account. In 
refs. [|3ll , [34|]3E|]3§fl the model equations can be valid for an OPO that is non degenerate either in polarisation or in 
frequency; however the intrinsic physical difference between the two cases was not made evident. In refs. p6| and ]5l| l 
the polarisation non-degenerate case and the polarisation and frequency non-dcgcncratc case have been considered 
and the role of the walk-off explored. Finally the case of polarisation coupling, due to cavity birefringence and/or 
dichroism has been analysed for the case of an OPO, showing Bloch walls formation |52||; these structures were later 
found also for second harmonic generation [p3[ . 

In this paper we will focus our attention on the polarisation pattern and vectorial defect formation in a Type II, 
frequency degenerate OPO. In particular we will consider the case in which the pump is not resonated in the cavity and 
can be eliminated from the dynamics by means of a multiple scale approximation. This simplification has also been 
made in order to take into account the effects of phase mismatch, which has been only partly studied in connection 
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with Type I OPO's |33]. In fact, the conservation of the photon momentum need not to be exact for down conversion 
to occur, the conversion efficiency being lower if k p — k s — hi 7^ 0. 

The paper is organized as follows. In section |n| we present the equations which model the Type II OPO and some 
exact, stable solutions, which are combinations of traveling waves (TW's). Due to the vectorial nature of the field these 
solutions represent spatial polarisation patterns. We also show that these TW's are spontaneously formed starting 
from a random perturbation of the trivial steady-state, that is represented by the field below the threshold of signal 



and idler generation. In section III we show that, on a background made of TW's, vectorial defects are spontaneously 
generated. The defects are classified as vectorial dislocations by comparing the quantities characteristic of a dislocation 
of a complex field with those of the defects obtained in the dynamical system. The dynamics of defects is also studied 



in this section. The stabilization of a single defect is then presented in section IV. The discussion of the effects 



due to the presence of the phase mismatch is given in section ^ and the conclusions presented in section VI. In 
Appendix I we derive the basic equations in which our study is based, and in Appendix II we derive the corresponding 
amplitude equation (a Swift-Hohenberg equation) for the instability leading to frequency-converted light generation 
in our model. 

II. EQUATIONS AND BACKGROUND SOLUTIONS 

We considered a model for a Type II OPO, for which the pump is not resonated in the cavity and in which some 
mismatch among the interacting waves is allowed. Previous studies of pattern formation, except ref. (33), neglected 
mismatch; however real devices are likely to present a residual mismatch due to the selection mechanism of the signal 
and idler frequencies (see [Q for details). Since the Type II configuration means that the signal and the idler are 
orthogonally polarized fields we will use the notation A x , A y to denote them. The starting point is the slowly varying 
envelope approximation propagation equations for a Type II interaction which include the mismatch (see for example 
p5[). By following the procedure outlined in Appendix I, we obtain the evolution equations for the cavity fields: 

d t A x = 7 X [-(1 + iA)A x + [iA* y + a\A y \ 2 A x ] + ia x V 2 A x (1) 
d t A y = 7 y [-(l + iA)A y + nA* + a\A x \ 2 A y ] + ia y V 2 A y (2) 

where j x -,"/ y are the decay rates of the signal and idler in the cavity, A the cavity detuning, a x ,a y the diffraction 
coefficients. Note that the detunings of the two fields are made equal to a common value A |5(| by choosing the 
temporal reference frame such that stationary homogeneous states do not have free rotating phases. The other terms 



which appear in these coupled equations are respectively: 

a) a linear parametric coupling with coefficient 

AAkL . sin(AkL/2) , . 

M = W — ) l kL/ i } (3) 

where Eq is the injected pump field, Ak the photon momentum mismatch and L the crystal length p7[ ; 

b) a cubic nonlinear coupling with coefficient 

a = 2 [(cos(AkL) - 1) + i(sin(AkL) - AkL)] (4) 

ZA/C Li 

At perfect phase-matching /i = Eq, cr = — 1 ; otherwise the coefficients are complex numbers; their real, imaginary 
parts are shown in fig. [l] for Eq — 1.5. Similar coefficients were also found for a Type I, singly resonant frequency 
degenerate OPO, by Longhi in ref. [[33]], with the same type of analysis. Note that the lasing of an OPO requires 
|Afc| < 7r/L because the parametric gain outside this range is too small and the device tends to "jump" to another pair 
of signal-idler frequencies which satisfy this condition |54j . Finally note that all coefficients can become dimensionless 
as soon as appropriate time, space and field amplitude scaling are performed. 

Equations (0,0) have exact traveling wave solutions of the same type as those found in ref. 63] : 



(a \ I 

X = v^ e iM+kT) 



(5) 



where 4>o is an arbitrary phase, while the intensity C and the phase difference a between the signal and the idler 
are given by 
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a 1 + a" A ± VImPM 2 ~ (c" - tr'A] 



(6) 



cosa = - a'C) + //'(A - a"C)\ (7) 

where a' = Re(cr),a" — Im(a), /i' = Re(fx),fi" = Im(fx) and A = A + |k| 2 (7 :E a 2 ; +J y a, y )/('Jx + ly) is an effective 
detuning parameter. The TW intensity C and the phase difference a are shown in fig. ^ for Eq = 1.5. 
The frequency shift oj is given by: 

^ = ^xMK- a ,)|k| 2 ] (8) 

This quantity can be zero if the diffraction coefficients are equal for the signal and the idler. Hereafter we will 
consider only the case of frequency degenerate OPO, i.e. when the frequency of the signal and idler are the same; this 
also implies that 7 X ~ "f y (without loss of generality we set ^ x — ^ y — 1 in our simulations, thus scaling time with the 



cavity lifetime) . The condition a x — a y — a can be exactly obtained for example by introducing compensating prisms 
in the cavity, as suggested in Appendix I. Strictly speaking if to = the solutions (|^) are not traveling waves but 
rather standing phase waves; however, we will generally refer to them as TW. It is important to note that considering 
the case u = is not restrictive in the sense that all features observed for this case survive also for u> ^ 0, at least 
when a uniform pump and periodic boundary conditions are used. We checked this by integrating eqs. also for 

a x =/= a y ; the main difference is that phase waves travel. Under the hypothesis made in this paragraph the effective 
detuning simply reads A = A + a|k| 2 . 

TW solutions (||) exist for |/i| 2 > 1 + A 2 = |/i c | 2 if a' + a" A < 0; in this case only the solution with the plus sign is 
acceptable in eq. (H). If a' + a" A > solutions exist for |//| 2 > |// c | 2 — ]7p( cr ' + c"A) and there is bistability up to /i c 
among the solutions obtained by taking the plus and minus signs. However, when no mismatch is present (Afc = 0), 
given that a' < 0, no bistability regime exists and TW are found only for \/^\ 2 > \/j, c \ 2 = 1. Note that if no mismatch 
is present the solutions also satisfy a = and therefore A x = A*. Hereafter we will consider only the case of perfect 
phase-matching and leave the observations about the effects of mismatch to section [y| 

A physical interpretation of the exact solutions can be given in terms of spatial polarisation patterns. In fact, under 
the assumptions made, the Stokes parameters |B8| 

s = \A x \ 2 + \A y \ 2 

5 1 = \A x \ 2 -\A y \ 2 

52 = Ax Ay + A x Ay 

5 3 = i(A x A* y - A* x A y ) (9) 
associated with the TWs at frequency u> g = u>i = w p /2 can be simply calculated as 

5 = 2C, Si = 0, S 2 = 2 C cos(2k c • r - a), S 3 = 2 C sin(-2k c • r + a) (10) 

These equations mean that the coupled- TW solution {A x , A y ) represents an elliptically polarized field with a two- 
valued spatially periodic azimuth (6 = tan -1 (S2 / Si) /2) and with a spatially periodic ellipticity parameter i](r) = 
a/2 — k c • r. In particular, in the direction of the wave- vector k c the polarisation changes from linear (with 9 = 7r/4), 
to right circular, to linear (with 9 = — tt/4), then to left circular and back again to linear (with 9 = ir/4) with a spatial 
period A c /2 = n/k c (see fig. ||). 

Other stable solutions can be found as combinations of the basic TW solutions (g) as shown in ref. |34| for a triply 
resonant, non-degenerate OPO. They can be generally written as: 
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A x \ 1 

= (hexp(ik c x) + f 2 exp(ik c y) + f 3 exp(-ik c x) + faexp(-ik c y)) | (11) 

The coefficients determine the kind of solution; for example the TW is (/i, /2, /3, A) = {VC, 0,0,0) and the 
standing rolls are (y/C, 0, s/C, 0). An interesting solution for this case is that of the alternating rolls (/i, /a, /3, Ja) = 
{VC,iVC, VC,iVC) i.e.: 

/ 



= 2VC{cos(k c x) + i cos(k c y)) (12) 



a; I V 1 



that gives rise to a square intensity pattern. We found that this structure is stable with periodic boundary conditions, 
but it deforms and disappears as soon as other, more realistic, boundary conditions are applied. We show in fig. ^ 
the square array with periodic boundary conditions after an evolution of t = 500 time units; no appreciable change 
of the initial condition ( |l2| ) has been observed. 

The equations also have a trivial steady-state solution A x — A y = 0, that represents the non lasing state of the OPO 
and whose stability analysis is similar to that performed in refs. |30|j3l|j35[] . Briefly, A x = A y = becomes unstable 
for \p\ > 1 if A < and for p > yl + A 2 if A > 0. For A > the trivial solution becomes unstable for homogeneous 
perturbations (k c = 0) and no pattern or structure is expected to be spontaneously formed. In numerical simulation 
we observed the transition to a final, stable, homogeneous state of intensity C, preceded by a transient regime where 
domain walls form but soon disappear jL7| . For A < the most unstable modes are TW with a critical wave- vector 
|k c | 2 = —A/a, i.e. such that the equivalent detuning A is zero. Note that the most unstable modes are also exact 
TW solutions of eqs. (|l]j|). Although stability of such TW's was not studied, numerical solutions give evidence that 
they are stable and moreover that they are spontaneously formed starting from a random perturbation of the trivial 
steady-state when the pump amplitude is above the threshold of the parametric interaction. Then, the spontaneously 



selected TW's have wave- vector modulus |k c | = k c — y/—A/a but a random orientation or phase in different region 
of space; this causes the occurrence of defects in the patterns as will be shown below. Hereafter we will treat only the 
case of polarization pattern formation, i.e. A < 0. 

If the pump amplitude is weakly above the threshold (5%) two different final situations can be found starting from 
a random perturbation of the trivial steady-state under periodic boundary conditions. In the first, all defects initially 
formed tend to annihilate each other and the final state is a single TW of random orientation, whose wave-vector is 
exactly k c . As in the case of non-degenerate OPO's |35| no intensity patterns appear but only phase stripe patterns, 
which correspond to the polarisation patterns described by eqs. (^|) and (|l^). In another case, for the same pump 



amplitude but with other initial random conditions, ordered structures of defects may form as shown in fig. o. By 
looking at the real and imaginary parts of the field A x it is clear that these defects are found along the fronts which 
separate regions where TW's with the same wave-vector but with different phase have been selected. From the figure 
it is also clear that these defects appear on a background of TW solutions, at those points where the stripes do not 
match. 

Defects also appear spontaneously with pump amplitudes well above threshold (50%); in this case the positions 
of the defects are not ordered, as shown for example in fig. ^ a) and b). Defects always appear in both orthogonal 
polarisation components (i.e. the signal and the idler simultaneously) and therefore they can be classified as vectorial, 
since the two components of the vector field (A x ,A y ) are zero at the same point in space. In fig. ^| c) and d) the 
phase of the polarisation components is also presented; note that it is not defined at the points where the amplitude 
goes to zero. Then, such objects require a topological classification that is furnished in the next section. Finally note 
that setting a x 7^ a y or y x 7^ ^ y does not alter the basic features observed; defects persist and are advected by the 
phase waves, which now travel as previously described. This has been cheked by running the same simulation (i.e. 
that of fig. ||) with and without equal coefficients. 

It is always useful to have in mind the general characteristics of the instability of the zero state in our OPO 
model, in order to relate it to other pattern forming systems. In our case, the instability for A < is at a non 
vanishing wavenumber fc c , so that a Swift-Hohcnberg equation for the complex envelope of the unstable modes is the 
appropriate amplitude-equation description close to threshold. For the case a x = a y , this Swift-Hohenberg equation 
has real coefficients. This equation is derived from our basic equations (1,2) in Appendix II. 

III. DEFECT CLASSIFICATION, FORMATION AND DYNAMICS 

In this section we analyse the formation and dynamics of defects, in Type II OPO's. The existence of defects in 
a Type I non-degenerate triply resonant OPO was pointed out in ref. pgj ] where the dynamics of an advected defect 
pair was also reported. Here, we present a detailed study of the defects which spontaneously appear in the Type II 
frequency degenerate OPO. The analysis includes: a) the classification of the defect type, showing that these defects 
are dislocations of the TW pattern; b) the formation process under several different conditions and regimes and c) 
the trapping and stabilization of a single defect. 

A dislocation of a scalar complex field A(r, 8), (r, 8) being spatial polar coordinates, can be defined as a configuration 
which can be deformed continuously (in a neibourhood of a point ro, which is then identified as the core of the 
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dislocation and it is taken as the origin of the polar coordinates) into the following function: 

A(r, 9) = R(r) exp(i(0 + kr cosO)) = R{r) exp(icf>) (13) 

The function R(r) is real and depends on the radial coordinate such that R(0) = 0. This condition is necessary 
because the phase <f> of the field is not defined in the origin of coordinates. In fig. ^ a) we show the real part of A from 
which we can observe that the defect is found at the point where the stripes of the background, whose wave-vector 
far from the defect is k, do not match. The phase <j> is shown in fig. |^ b); the defect is centered at the point ro = 0, 
around which the phase changes by 2tt. The topological charge associated with the phase change is defined by the 
integral of the phase gradient on a closed contour surrounding the point ro 

no = ^~ <f V<£ • dr = 1 (14) 

Note that there exists also another dislocation obtained by setting 9^—8 for which the topological charge no is 
equal to —1. However the most interesting feature of a dislocation is that it is characterised by another topological 
charge. This is associated with the function 

•"^ = ""•"' ({ft) (15) 

where the indices x, y of the term on the right hand side refer to the vector components. This function, which 
defines the orientation angle of the wave- vector k, is shown in fig. [?] c) for the field A defined by eq. (|l3|). Note 
that there are two points where ip is not defined: one is located at the point ro, corresponding to the point of zero 
amplitude; the second is at another point ri where \r± — ro\ = k^ 1 . Then, the topological charges associated with ip 



too = — <P V-0 • dr = 1 (16) 

where the integration contour surrounds the point ro but leaves outside the point ri and 

mi = — / W • dr = -1 (17) 
2tt J ri 

where the integration contour surrounds the point ri but leaves outside the point ro . Note that there is not a zero 
in the amplitude R at the point ri . 

We now turn our attention back to the Type II OPO and compare the defects which are found in the numerical 
solutions of eqs. (0,0) with the dislocation just defined. In this section, to better compare numerical and theoretical 



results only periodic boundary conditions and homogeneous pump beams are considered; more realistic conditions 



will be explored in section (IV) when stabilization of a dislocation is considered. 

The randomly generated defects of figures || and || are dislocations on a background of TW's, although in the 
second case it is difficult to recognize the background solutions since the defects are dense in the integration window. 
However a fingerprint of TW is found also for the data of fig. |^ by observing the far field in fig. ||. There, a ring of 
modes of radius k c is clearly observable, a clear indication that the background is made mostly of TW solutions of 
different orientation but the same wave- number k c . The dislocations can be also clearly identified when the quantities 
tfj x and tfjy (defined again through eq. ( ^5|) where the phases are those of A x ,A y ) are calculated for the same data of 
fig. |(| The result is presented in figures || a) and b) where a contour plot of \A X \, \ A y \ taken at a low intensity level 
to indicate the locations of the zeros of the field, is superimposed on the grey-scale coded images of the functions ip x 
and ip y . From fig. ^| we note that in this static condition "chains" of defects of the functions ip x ,y are formed. Two 
zeros of the amplitude (defects of charges no, mo) are in fact connected by the relative pair of defects of the functions 
if> XtV at the second point (defects of charge mi). This stems from the fact that the phase difference between A x and 
A y is zero (when there is no mismatch) and thus A x = A*. 

Regarding the dynamics |5{| of the spontaneously generated defects, we observe that they appear at the early stage 
of the process. Some of them are annihilated within a few thousand cavity lifetime units. After this stage we cannot 
observe any annihilation, up to 100000 cavity lifetime units. The initial number of dislocations and their positions 
change with different initial conditions. The statistical distribution of the final number of defects, by changing initial 
conditions, is shown in fig. [H]; observe that a Gaussian fit is actually very good. We note that the defect density 
never reaches that of the square pattern (eq. (|l2|)) with critical wave- vector k c . After the initial formation stage, we 
observe in our runs that defects have a slow, random motion. This regime is mostly dominated by the attraction of 
pairs of defects of opposite charges rig = ±1. For longer times the motion tends to stop and an equilibrium condition 
is reached. Although no detailed statistical study has been conducted we also note, by inspection of fig. Q that there 
seems to exist a characteristic length of separation between the defects, which is an indication that an equilibrium is 
likely to be reached. 

To better study defect interaction we use ( |l3| ) as an initial condition and integrate eqs. (Q,0). Due to the periodic 



boundary conditions a pair of defects is generated, as shown in fig. 11a). These correspond to two dislocations in the 
real part of each field (fig. |ll]b)), two phase singularities with opposite topological charge (fig. |ll] c)) and two pairs 
of defects associated with the quantity ijj x (fig. |ll| d)). At the initial stage the zeros of the amplitude, which have 
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opposite no charge, attract each other; later, when the ipx defects (which have the same charge mo = 1) get closer, 
this attraction is stopped and an equilibrium position is found p9| . This paired structure is stable at least up to the 
time explored by numerical solutions. 

In general, the detailed behavior of defects dynamics is governed both by local, phase curvature effects and by 
defect-defect interactions. The range of these interactions depends on the parameters of the system, which determine 
the spatial extent of the defect. For example if we make the detuning more negative (e.g. A x — A y — — 1 and 
a x = CLy = 1) two defects displaced parallel to the background wavefronts either annihilate, if close enough, or repel 
each other, if further away. 

Figure |l2| also shows a pair of defects (same parameter values as in Figure |ll| but different initial displacement) 
undergoing motion; the background has a wave-vector of 2k c and it is stable when the system is pumped far above 
threshold. When their separation is small enough, the defects briefly rotate around each other before continuing their 
transverse motion. As a consequence of this motion and the periodic boundary conditions, the defects collide and 
annihilate shortly afterwards (not shown). 

Figure ^| shows another interesting example. Two defects approach each other while the background phase wave 
undergoes a zig-zag instability [ |60[ . When they meet they form a bound pair which follows the local curvature of the 
background phase fronts, changing direction as the slow phase dynamics alters this curvature. The bound defect pair 
is stable for the duration of the simulation (115000 cavity lifetimes). It also exists for lower pump powers, down to at 
least five percent above threshold. On the other hand, at slightly higher pump powers than that used in Figure [l^, 
the structure destabilises: after a short time the defect pair rotates slightly before annihilating. 

The previous examples show that long-lived defect interactions are important features of the dynamics of Type 
II OPO's. This fact is in agreement with instances of long term survival of many defects generated by random 
initial conditions (e.g. Fig. ^J) and is clearly different from vortex dynamics observed, for example, in complex 
Ginzburg-Landau models ]6l| ] . 

Close to threshold, it is possible to study the behavior of the Type II OPO dislocations by means of amplitude 
equations. These are derived in Appendix II where comparisons are made with simulations of the full model. 

IV. TRAPPING OF ISOLATED DISLOCATIONS 



We know that periodic boundary conditions always require the total topological charge within the domain to be 
zero. In this section, we take into account that the pump amplitude Eq actually has a spatial dependence and periodic 
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boundary conditions are not applicable. If a spatial Gaussian or super-Gaussian distribution is used for the parametric 
gain factor [i (eq. (0)), the spontaneously generated defects often tend to move to the region where the field is zero. 
The trapping of single defects necessary for their experimental observation as single entities is then an interesting 
issue. 

We tried, successfully, to isolate a single defect in a more realistic pump beam, by using, as an initial condition for 
one of the fields, a doughnut mode. This idea follows the experimental results of induced dislocation formation in 
quadratic nonlinear interactions [ |62| . Note that no cavity was used in those experiments and the defect was maintained 
by forcing a doughnut mode in the field at the crystal input. Here the dislocation is "written" at the beginning but, 
after that, no injection is required since the cavity provides the necessary feedback to sustain the structure. The 
resulting, stable trapped defect is shown in fig. |l4|. The spiral wave in the phase field spins when a x ^ a y , otherwise 
it is stationary after the initial formation stage. Note that in order to trap the defect the beam size must be kept 
smaller than the wavelength of the most unstable TW A c = 2ir/k c otherwise more defects are induced or the initial 
defect is dragged outside the beam. The creation of other defects is due to the tendency of the system to generate 
pairs. If the system size is comparable with the critical wavelength, the paired defect will appear inside the pumping 
region and will destabilize the original defect. If the size is smaller the second defect will be generated outside the 
pumping region and therefore will not influence the dynamics of the trapped defect (fig[l4| c)). The amplitude profiles 
of the fields of this single defect are shown in fig. [15[ 

This way of generating the defects is also useful for understanding the mechanism of their formation. The presence 
of a defect in one field, say A x , induces a defect of opposite charge no in the other because of the parametric coupling. 
In fact if A x ~ at t = the term which dominates the dynamics of A y in the initial regime is the largest linear one, 
i.e. [iA* x . Hence, A y is forced by the complex conjugate of A x and a defect of opposite charge is formed. Eventually the 
other terms of eqs. (fj]j|) become significant: the real part of the cubic term provides the saturation which stabilizes 
the solution and nonlinear phase modulations appear if mismatch is present, due to its imaginary part. 

If we initialise a defect with topological charge no = ±2 using the technique described above, we observe that it 
breaks up into two defects of charge n = ±1. 

Since the generation of such defects has been seeded externally, the question arises whether it is possible to obtain 
defects for positive detunings. In spite of the absence of stable traveling-wave solutions, single vortices can be trapped 
by using the same seeding technique. 

We also observe that walk-off removes all structures from the pumping region and a single TW is selected asymp- 
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totically. However, it is interesting to note that in the transient regime dislocations are formed. They appear at the 
front which divides two regions where stripes with different wave-vector modulus are selected. The difference in the 
wave- vectors is due to the walk-off, as shown in ref. jl6| . 

Finally, we estimate the threshold for the observation of this isolated dislocation in an OPO. Considering data 
for KTiOPOi [Q: \ — 7.33pm/V,no = 1.8 at Ao = 1.064/xm, and by using the definition of the effective cou- 
pling paramater defined in Appendix 1 (a) and considering a nonlinear crystal length of 1cm, a mirror transmit- 
tivity 1CP 2 , the input field at threshold (E=l) can be evaluated as Ejn,o — 2l0 4 V/m that yields an intensity 
I = no/2 (/Lto/eo) 1 ' 2 |£'/iv,o| 2 — lMW/m 2 (fio, eq are respectively the vacuum magnetic permeability and electric per- 
mittivity constants). The super-Gaussian beam used in the simulation has a beam diameter of about 10 normalized 
units; a spatial normalized unit corresponds, for the diffraction coefficient of the simulation to about 0.195mm; hence 
the CW threshold power is about 0.95VF. The physical cavity decay rates -y x y are about Q.HGHz while the cavity 
detuning is SAM Hz. All these data are compatible with an experimental realization. 



V. MISMATCH EFFECTS 

This section is dedicated to the study of the effects of mismatch (Ak ^ 0) between the fields in the OPO. As 
previously noted, in real devices it is often possible that the selection rules of the oscillation frequencies force the 
OPO to emit radiation with a slight mismatch |j54| . 

The effects of the mismatch, predicted by the analysis of section ||, are the following: a) an increase in the threshold 
for signal generation (i.e. the instability, see eq. (|!|)); b) a spatial shift of the iso-polarisation lines for the exact TW 
solutions. The latter effect is due to the appearance of the phase shift a ^ (see eq. (|^)) between polarisation 
components; this stems from the imaginary parts of the coefficients fi and a, the contribution of the latter being a 
nonlinear phase modulation due to the x' 2 ' : cascading effect [Q. 

For the case of the spontaneous generation of defects we have been able to determine that the phase shift a is 
non zero in the regions where the background solution dominates. The numerically found value is very close to the 
value predicted by eq. (0). However, approaching the zeros of the amplitude the phase difference also tends to zero 
although the phase of the fields is not strictly defined at the defects. In practice the zeros of the spatial function 
a(x, y) are located exactly at the positions of the field defects. This behavior can be explained as follows: close to 
the defects the cubic terms tend to zero; in particular the term proportional to the imaginary part of a (i.e. the 
Kerr-like phase modulation term) is the first that can be neglected. Although /i = \/i\exp(i/3) is complex, by an 
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appropriate re-definition of the variables as Aj = Ajexp(i(3/2) j = 1, 2 we can obtain a set of equations for Aj where 
the parametric gain is purely real (//' = 0). Thus, we reduce the problem to the case where no mismatch is present 
(see eq. (^)) and therefore a — > 0. 

The final number of spontaneously generated dislocations seems also to be affected by the phase mismatch. We 
note that this number decreases by changing the mismatch from negative to positive values, keeping all the other 
parameters fixed. In particular it seems that positive values of the mismatch encourage the annihilation of the pairs 
of dislocations. This has been verified by numerical integrations with the same parameters and initial conditions as 
fig. [ll] but with non zero mismatch. For Afc > we observed that the pair is annihilated, while for Afc < the 
pair reaches a stable configuration. A possible cause of this asymmetry is the phase modulation a introduced by the 
mismatch, which changes sign according to the sign of the mismatch (see fig. |J). 

The effect of phase mismatch on the trapped defect is to destabilise it, by generating asymmetries in the fields for 
both positive and negative values of Afc. 



VI. CONCLUSIONS 

In conclusion we have studied the polarisation pattern and vectorial defect formation in a Type II frequency 
degenerate optical parametric oscillator. 

We found that the preferred solutions, those that are selected out of an initial perturbation of the zero state, are 
conjugate traveling waves in the two components of the polarisation. A spatial polarisation pattern is then formed: 
the state of polarisation changes along the spatial coordinate parallel to the selected wave- vector k c with period equal 
to A c /2 = ir/k c . In particular the state of polarisation changes along the meridian of the Poincare sphere which 
passes through the linearly polarized states of azimuth 8 = ±ir/4 if Afc = 0. The magnitude of the wave-vector of 
the selected solution is fixed by the parameters and has a single random orientation close to threshold and multiple 
random orientations far from threshold. Combinations of traveling waves, that form square patterns are also stable 
solutions and, close to the threshold of the instability, ordered arrays of defects can form spontaneously. 

Such defects are isolated zeros of the two linear components of the polarisation, i.e. they are vectorial defects. 
They are dislocations which form in spatial positions where the background solutions (traveling waves) do not match 
spatially. Two different kinds of topological charge must be defined: one kind of charge is associated with the phase 
and two charges with the director angle of the field wave-vector. The first charge, located at the point where the 
amplitude goes to zero, can be ±1 and is always opposite in the two polarisation components. The charge associated 
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with the wave- vector is always +1 at the point where the amplitude is zero and —1 at the paired point, where the 
amplitude is not zero. The polarisation components are totally correlated; in fact all defects, both in the phase and 
wave-vector fields, have a corresponding defect in the other polarisation component. In this way the defects form 
chains in which the separation among defects seems to have a typical size, which is of the order of the background 
wavelength. 

The trapping of an isolated defect has also been demonstrated. This is accomplished by keeping the size of the 
pump beam smaller than the critical wavelength of the preferred traveling wave such that a second defect cannot be 
created inside the pump beam but rather lies outside and does not influence the dynamics. 

Finally, we have addressed the effects which may arise when the nonlinear interaction is slightly phase mismatched, 
i.e. the increase in the threshold of signal generation, and the linear and nonlinear phase shifts among the polarisations, 
the latter due to the cascading effect. Numerical results showed that a positive mismatch can favor the annihilation 
of dislocations. 



Possible applications of defects can be foreseen in the field of particle and atom trapping 65 
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VII. APPENDIX I 

In this appendix we derive the dynamical equations for a cw Type-II OPO in the presence of diffraction. We start 
from the amplitude equations for the pump, signal and idler fields in the crystal [p5fl 



— d t E = —V 2 E 5£ Ex E y e- lAkz 

C ZKq TIqC 



d z E x + '± d t E x = ^ 2 E X + E E* (18) 



^ dtEx = ^-^E x + ^x e e; 

d z E y + %L dt E y = ^ 2 E y + E Q E* e iAk *. 



where _E , E x , E y are the slowly varying amplitudes of pump, signal and idler respectively, J7 = £l x + £l y is the 
frequency constraint on the OPO, fcj = nfili/c are the wave-numbers, Afc ~ kg — k x — k y is the phase mismatch, 
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X is the second order susceptibility of the crystal, and c is the speed of light. We treat here the perfectly matched 
case (Afc = 0) in order to focus on the form of the diffraction coefficients of the final equations. The phase matching 
condition implies that once the three frequencies and two refractive indices are given, then the third refractive index 
is determined. For example for no we have 

^X^X ~t~ TlyQy Tl X -(- Tiy (l9) 

^^x ^^"y 2 

where the second equation is valid at frequency degeneracy. It is also useful to introduce Q x = /if2o, Q y = vSl® with 
H + v = 1 and the effective coupling parameter a — 47rf2ox/( nc ) to obtain 



d z E Q + — d t Eo = — V 2 £ - oiE x Ey 



c 2kq 

d z E x + ^ d t E x = 7^r^ 2 Ex + liaEoE* y (20) 



d z E v + -s- d t E y = —V 2 E y + vaEoE*. 



All these equations describing the fields in the crystal have a similar form of the kind 

d z Ei + —d t E t = idiV 2 Ei + & NLT(Ej, E m ) (21) 

where NLT stands for Non Linear Terms and the indices i, j, m take the symbolic values 0, x, y. 

We consider a ring cavity of legth C with a crystal length L and apply longitudinal boundary conditions and 
the standard Mean Field Limit (MFL). We consider here the triply resonant case, the doubly-resonant case of non 
resonated pump being a subset of this general frame. We then obtain 

+ r I 7— TTr d * F * = ~ ^ + *OiV 2 £< + liEiN,i + PiNLT(Ej, E m ) (22) 

£ + (n, - 1)L 

where 



t' = t + 



C-L 



c 



z 
L 



coi e - cSi (oji - Qi)£ 

7t = "7> I 7 TTT! d * ~ 



CC ft= ,, # ,, (23) 



C+(m-i)V £ + (m-i)L c+(m-i)V 

c£ , a _ c(3. t L 

ai ~ 2ki[£ + (m - 1)L] ' Pi ~ L + (m - l)L' 

Note that the presence of birefringence introduces an explicit dependence on the refractive indices in the coefficients 
of the final equations. This is a consequence of the tight comb of resonances observed when the length of the cavity 
is scanned |56| . 

Finally, we introduce a linear transformation of the fields 
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and a final normalisation of the parameters 



The final equations read 



E(x,y) = — J E INfi Ai = — = . (25) 



ft' A) = 7o [-(1 + «A )A + E(x, y) - A x A y ] + ia Q V 2 A a 

d t >A x = j x [-(1 + iA x )A x + A Al] + ia x V 2 A x (26) 
dfAy = [-(1 + iA y )A y + A A%] + ia y V 2 A y 

where 

= g £ = c£ 

7l £ + (7i 4 - ° l 2fci[/: + (n< - 1)1] ' 1 ' 

These equations clearly show that the loss and diffraction coefficients depend critically on the refractive indices even 

in the frequency degenerate case. 

Since our analysis starts from stationary homogeneous solutions, the choice of the temporal reference frame is fixed 

by the condition = A y = A which excludes phase rotations for the stationary homogeneous states p6[ . In this 

case the final equations are 

dt'Ao = 70 [-(1 + iA Q )A a + E(x, y) - A x A y ] + ia V 2 A 

d t >A x = lx [-(1 + iA)A x + A A*] + ia x V 2 A x (28) 
d t .Ay = ly [-(1 + iA)A y + A A*} + ia y V 2 A y . 

We note that it is still possible to have equal loss and diffraction coefficients for the three waves if we consider 
cavities of different lengths Ci for each field. This can be achieved by inserting compensating prisms of chosen length 
Ai and refractive indices iVj. In this case the equations are the same as ( p8| ) but with redefined loss and diffraction 
coefficients 

~. = °-^i± • a = °£ (29) 

h d + (m - l)L + (Nt - 1) A, ' 1 2k i [C i + (n i -l)L+(N i -l)A i ]' v ' 

By adjusting losses <Xj and compensating prism coefficients, one can select equal 7^. In this way the compensating 

prism parameters can be left free to adjust the diffraction coefficients dj. In particular, we want a x = a y at degeneracy 

which means 
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n x [C x + (n x - l)L + {N x - 1) A x ] = n y [C y + (n y - l)L + {N y - 1) A y ] (30) 

which leaves a lot of flexibility for the final setting. 

Finally, we observe that in the absence of a cavity for the pump field there are just two equations for the signal 
and idler waves remaining. By setting Aq = E — A x A y in equations (Pq) , we obtain equations for the oppositely 
polarized fields A x and A y in a doubly-resonant configuration. 



VIII. APPENDIX II 

In this Appendix we perform a weakly nonlinear analysis in the region k c <C 1 and close to threshold in order to 
determine the appropriate amplitude-equation description for the instability. Related derivations have been presented 
in |3q»3] • For simplicity, consider the case 



J x = 7 B = 7 



(31) 



With the benefit of some a posteriori knowledge, we rewrite the field equations in terms of the variables F\ — A x + A* 
and F 2 = A x - A*: 



d t Fi = 7 
d t F 2 = 7 



1 



(E-1)F 1 --F*{F 



-{E+l)F 2 + -FZ(F^Fl) 



+ i (aV 2 - 7 A) F 2 
+ i(aV 2 - 7 A)Fi. 



We define a smallness parameter e and perform appropriate scalings and expansions: 



(32) 



E = 1 + e 2 P , T — e t , X = sfe x , Y = y/e y , A = eS , F l = eF a + e 2 F i2 + 



(33) 



We can then proceed order by order in e in a straightforward manner. At third order we have 

2 



d t F x = 7 



(E-1)F 1 -~\F 1 \ 2 F 1 



F 2 = irv 2 -A)r, 



(34) 
(35) 



The important point is that there is a single order parameter (F\) governed by a Swift- Hohenberg equation (SHE) 
fiof and a second field (F 2 ) which is slaved to F\. In other words, the dynamics of the system is described by a 
single complex scalar field. Note that it is only when A x = A* that the order parameter becomes real |3£|. An extra 
imaginary transverse Laplacian term appears in (52) when a x ^ a y . 
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Despite the assumptions underlying its derivation (Equations (p3|)), the SHE can sometimes offer a good description 
of the system even outside its expected range of validity. For example when integrating eq. ( |34| ) with the same 
parameter values as in Figure O, the final state is essentially the same stable defect pair. Nonetheless, notable 
failures of the SHE can be found; for example in reproducing the movement of the bound defect pair in Figure [D| 
This indicates that the original dynamics is not well reproduced when the assumptions leading to ( |33| , |34| ) are violated, 
and serves as a warning about the limitations of order parameter equations such as . 



+ http:/ /www. imedea.uib.es/PhysDept 

[1] Nonlinear optical structures, patterns, chaos, ed. L.A. Lugiato special issue of Chaos, solitons and fractals 4, 1251 (1994) 
and references therein. 

[2] Nonlinear optical structures, patterns, chaos, ed. L.A. Lugiato special issue of Chaos, solitons and fractals 10, 617 (1999) 
and references therein. 

[3] Proceedings of the Euroconference: Patterns in nonlinear optical systems, Alicante (Spain) 1998 published in the special 
issue of Journal of Optics B: Quantum Semiclass. Opt., 1, (1999). 

[4] M. Tlidi, P. Mandel, R. Lefever, Phys. Rev. Lett. 73, 640 (1994). 

[5] M. Saffman, D. Montgomery, D. Z. Anderson, Opt. Lett. 19, 518 (1994). 

[6] W. Firth, A.J. Scroggie, Phys. Rev. Lett. 76, 1623 (1996). 

[7] N. N. Rosanov, Progress in Optics 35, 1 (1996). 

[8] A. Schreiber, B. Thuring, M. Kreuzer, T. Tschudi, Opt. Comm. 136, 415 (1997). 

[9] C. Etrich, U. Peschel, F. Lederer, Phys. Rev. Lett. 79, 2454 (1997). 

[10] K. Staliunas, J. V. Sanchez-Morcillo, Phys. Rev. A 57, 1454 (1998). 

[11] L. Spinelli, G. Tissoni, M. Brambilla, F. Prati, L. Lugiato, Phys. Rev. A 58, 2542 (1998). 

[12] C. O. Weiss, M. Vaupel, K. Staliunas, G. Slekys, V. B. Taranenko, Appl. Phys. B 68, 151 (1999). 

[13] G.-L. Oppo, A.J. Scroggie, W.J. Firth, Journal of Optics B: Quantum Semiclass. Opt. 1, 133 (1999). 

18 



[14] M. Le Berre, D. Leduc, E. Ressayre, A. Tallet, Journal of Optics B: Quantum Semiclass. Opt., 1, 153 (1999). 
[15] B. Schapers, M. Feldmann, T. Ackemann, W. Lange, Phys. Rev. Lett, 85, 748 (2000). 

[16] P.L. Ramazza, E. Benkler, S. Boccaletti, U. Bortolozzo, S. Ducci, F.T. Arecchi, "Tuning the properties of localized 



structures in a nonlinear interferometer" (Los Alamos e-print nlin.PS/0109016). 
[17] G.L. Oppo, A.J. Scroggie, W.J. Firth, Phys. Rev. E, 63, 066209 (2001). 
[18] J.B. Geddes, J.V. Moloney, E.M. Wright, W.J. Firth Opt. Comm., Ill, 623 (1994). 
[19] M. Hoyuelos, P. Colet, M. San Miguel, D. Walgraef Phys. Rev. £58, 2992 (1998). 
[20] R. Gallego, M. San Miguel, R. Toral, Phys. Rev. £61, 2241 (2000). 

[21] A. Aumann, E. Bthe, Yu. A. Logvin, T. Ackemann, W. Lange, Phys. Rev. A, 56, R1709 (1997). 
[22] M. San Miguel, Phys. Rev. Lett. 75, 425 (1995). 

[23] J.Martin-Regalado, S. Balle, M. San Miguel, Optics Letters 22, 460 (1997). 

[24] E. Hernandez-Garci'a, M. Hoyuelos, P. Colet, R. Montagne, M. San Miguel, Int. J. Bif. and Chaos 9, 2257 (1999); E. 
Hernandez-Garci'a, M. Hoyuelos, P. Colet, M. San Miguel, Physical Review Letters, 85, 744 (2000). 

[25] U. Peschel, D. Michaelis, C. Etrich, F. Lederer Phys. Rev. E 58, R2745 (1998). 

[26] P. Coullet, L. Gil, F. Rocca, Opt. Comm., 73, 403 (1989). 

[27] A.B. Coates et al., Phys. Rev. A, 49, 1452 (1994). 

[28] F.T. Arecchi, G. Giacomelli, P. L. Ramazza, S. Residori, Phys. Rev. Lett., 67, 3751 (1991). 
[29] K. Staliunas, G.Slekys, C. O. Weiss, Phys. Rev. Lett, 79, 2658 (1997). 
[30] G-L. Oppo, M. Brambilla, L.A. Lugiato, Phys. Rev. A, 49, 2028 (1994). 

[31] G-L. Oppo, M. Brambilla, D. Camesasca, A. Gatti, L.A. Lugiato, J. Mod. Opt. 41, 1151 (1994). 
[32] K. Staliunas, J. Mod. Opt. 42, 1261 (1995). 
[33] S. Longhi, J. Mod. Opt. 43, 1089 (1996). 
[34] S. Longhi, J. Mod. Opt. 43, 1569 (1996). 

19 



[35] S. Longhi, Phys. Rev. A 53, 4488 (1996). 

[36] S. Longhi, A. Geraci Phys. Rev. A 54, 4581 (1996). 

[37] G.J. de Varcarcel, K. Staliunas, E. Roldan, V.J. Sanchez-Morcillo, Phys. Rev. A 54, 1609 (1996). 

[38] V.J. Sanchez-Morcillo, E. Roldan, G.J. de Varcarcel, K. Staliunas, Phys. Rev. A 56, 3237 (1997). 

[39] M. Santagiustina, P. Colet, M. San Miguel, D. Walgrael, Phys. Rev. E 58, 3843-3853 (1998). 

[40] H. Ward, M.N. Ouarzazi, M. Taki, P. Glorieux, Eur. Phys. Journ. D. 3, 275-288 (1998). 

[41] P. Lodahl, M. Bache, M. Saffman, Phys Rev. Lett. 85, 4506 (2000). 

[42] S. Trillo, M. Haelterman, A. Sheppard, Opt. Lett., 22, 1454-1456 (1997). 

[43] S. Longhi, Phys. Scr. 57, 611 (1997). 

[44] S. Trillo, M. Haelterman, Opt. Lett. 23, 1514-1526 (1998). 

[45] M. Santagiustina, P. Colet, M. San Miguel, D. Walgrael, Opt. Lett. 23, 1167-1169 (1998). 

[46] G. Izus, M. Santagiustina, M. San Miguel, P. Colet, Jour. Opt. Soc. Am. 16, 1592 (1999). 

[47] M. Vaupel, A. Matre, C. Fabre, Phys. Rev. Lett., 83, 5278 (1999). 

[48] S. Ducci, N. Treps, A. Maitre, C. Fabre, Phys. Rev. A, 64, 023803 (2001). 

[49] M. Santagiustina, P. Colet, M. San Miguel, D. Walgrael, Opt. Express 3, 63-70 (1998). 

[50] M. Taki,N. Ouarzazi, H. Ward, P. Glorieux, Jour. Opt. Soc. Am. B, 17, 997 (2000). 

[51] H. Ward, M. N. Ouarzazi, M. Taki, P. Glorieux, Phys. Rev. A, 63, 016604 (2001). 

[52] G. Izus, M. Santagiustina, M. San Miguel, Opt. Lett, 25, 1454 (2000); G. Izus, M. Santagiustina, M. San Miguel "Phase- 
Locked Spatial Domains and Bloch Domain Walls in Type-II Optical Parametric Oscillators", to be published, Physical 
Review E (2001). 

[53] D. Michaelis, U. Peschel, F. Lederer, D. V. Skryabin, W. J. Firth, Phys. Rev. E, 63, 066602 (2001). 
[54] R.C. Eckart, CD. Nabors, W.J. Kozlovsky, R.L. Byer, Jour. Opt. Soc. Am. B 8, 646 (1991). 

[55] J-Y Zhang, J. Huang, Y.R. Shen, "Optical parametric generation and amplification", in Laser Science and Technology vol 

20 



19, V.S. Letokhov, C.V. Shank, Y.R. Shen, H. Walther Eds. (Harwood Ac, 1995). 
[56] T. Debuisschert, A. Sizmann, E. Giacobino, C. Fabre, J. Opt. Soc. Am. B 10, 1668 (1993). 

[57] In our approximation the cavity is almost filled by the nonlinear crystal and thus the cavity length is very close to the 
crystal length. 

[58] B. E. A. Saleh, M. C. Teich, "The fundamentals of photonics" , (Wil ey, New York 1991), Chap. 6. 

[59] Movies showing time evolution under several circumstances can be downloaded from the URL 
http : //www. imedea.uib . es/Nonlinear/research_topics/OPO/. 

[60] M.C. Cross, RC. Hohenberg, Reviews of Modern Physics 65, 851 (1993). 

[61] L.M. Pismen, Vortices in Nonlinear Fields, (Oxford University Press, Oxford, 1999). 

[62] A. Berzanskis, A. Matijosius, A. Piskarkas, V. Smilgevicius, A. Stabinis,Opt Comm. 140, 273 (1997); Opt. Comm. 150, 
372 (1998). 

[63] C.L. Tang, L. K. Cheng, Laser Science and Technology vol 20, V.S. Letokhov, C.V. Shank, Y.R. Shen, H. Walther Eds. 
(Harwood Ac, 1995). 

[64] G.I. Stegeman, M. Sheik-Bahae, E. Van Stryland, G. Assanto, Opt. Lett. 18, 13 (1993). 
[65] K.T. Gahagan, G.A. Swartzlander Jr.,Op£. Lett. 21, 827 (1996). 



21 



FIGURE CAPTIONS: 

Fig. 0: The solid (dashed) curve is the real (imaginary) part of /x and the dashed-dotted (dotted) curve the real 
(imaginary) part of a as functions of the mismatch (Eq = 1.5). 

Fig. ||: The solid (dashed) curve is the traveling wave solution intensity C (phase difference a) as a function of the 
mismatch, as given by eqs. (||J7|) for A = 0. 

Fig. |j: The change of the state of polarisation of the solution (^) (k c = (fc c , 0), C — 0.8, a — 0, 4>q = 0) is shown at 
the spatial points of coordinates x ~ nir/8, n — 0..7 (left to right and top to bottom). 

Fig. |: The amplitude (\A X \) of the alternating roll solution (|||) at t = 500 time units (a = 0.8, A = -0.8, E = 1.5, 
the integrating window size was Xl =Yl = 62.8 units). 

Fig. |{5|: Ordered arrays of dislocations obtained close to threshold (E = 1.05): a) \A X \; b) Re(A x ); c) Im(A x ). The 
other parameters were a x = 2, a y = 2.05, ^ x = 1, j y = 1.025, A = -0.8, X L =Y L = 80, t = 4000. 

Fig. H Amplitude and phase of the polarisation components of the field well above threshold (Eq — 1.5, a = 1, A = 
-0.25, X L = Y L = 80, i = 2000). a) \A X \; b) \A y \: c) tan~ 1 [Im{A x ) / Re{A x )]] d) tan- 1 [Im{A v ) / Re(A v )]. Note that 
the defects are always in pairs and that signs of the charge uq are opposite in the two components. 

Fig. 0: Dislocation of a complex field: a) Re(A); b) cf>; c) ip; where A, <fi are defined by eq. (13) and tp by eq. ( |l5| ) 
(fc = l,i?(0) = 0,i?(oo) = 1). 

Fig. ||: The amplitude of the far field (Fourier transform) of A x for the same data of fig. |6[ The average modulus 



of the wavenumber of the modes in the ring is about 0.47 and the predicted value is k c = y— A/a = 0.5. 

Fig. [| Functions a) ip x and b) tp y for the data of fig. || Contour plots, with a contour level close to zero, of 
\A X \, \ A y \ have been superimposed to locate the positions of the amplitude zeros. 

Fig. [lo| Distribution (bars) of the final number of defects for different random initial conditions (a = 1,A = 
—0.4, Eq = 1.5, Xl = Yl = 50). The Gaussian fit average and variance are those obtained form the data. 

Fig. [ll]: Pair of dislocations obtained by integrating (0,||) with a single dislocation ( |l3| ) as an initial condition 
for A x and a random initial condition for A y : a) | ^4.^ | ; b) Re(A x ); c) 4> x ; d) ip x . The integration parameters were: 
a = 2, A = -0.2, E a = 1.5, X L = Y L = 100, t = 2000 

Fig. [l2]: Interaction of two vortices (labelled 1 and 2) for j x = 7^ = 1, A = —0.2 and a x = a y = 2, Xl = Yl = 
16ir/k c . The images show \Ai\ at intervals of 50 cavity lifetimes. 

Fig. |l3|: A bound defect pair for j x = *y y = 1, A = — 1 and a x = a y = 1, Xl = Yl = 32ir/k c . Figures (a), (c), (e) 
and (g) show the real part of A x at t=5500, 11000, 36000 and 54000 cavity lifetimes respectively. The figures in the 
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rightmost column show \A X \ at the corresponding times. The arrows indicate the direction of transverse motion of 
the defects. 

Fig. [li]: Isolated defect: a) \A X \; b) <ft x ; c) ip x . The initial condition for A x was a doughnut mode of the form 
A x (x,y) = 0.01 (x + iy) exp(—(x 2 + y 2 )/200) while A y was random. The parameters of the integration were: 
a = 2, A = -0.2, E = 1.5, X L = Y L = 60, t = 2000. 

Fig. 15: Amplitude profile of the isolated defect of fig. |l4|: \A X \ solid curve, \A y \ dashed curve. 
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